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In this report we propose that a charge rearrangement within heavy ions may occur during fusion 
reactions, lowering the Coulomb barrier as the nuclear systems approach each other. In this situation 
a modified Gamow factor has been calculated, which emphasizes the importance of spatial degrees 
of freedom in nuclear matter. 
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I. INTRODUCTION 



It is commonly believed that the point-like approxima- 
tion for the Coulomb potential is accurate enough when 
describing the nuclear dynamics during a fusion interac- 
tion In this context, there are a number of anoma- 
lies observed in fusion reactions near the Coulomb barrier 
and a number of anomalies in the abundances of 
p- elements produced in supernovas specially the pro- 
duction of rare nuclide, that do not fit well in our present 
understanding of supernova and nuclear dynamics. 

In this work, we will consider the effect of a spa- 
tial redistribution of charge during a nuclear interaction, 
which may add an extra ingredient when considering the 
anomalies observed in fusion reactions near the Coulomb 
barrier 011]. The key idea in this model is to take each 
nuclei as a plasma hall (inspired in a quark-gluon plasma 
although no calculation presented here has any re- 
lation to that specific subject), that has the ability to 
redistribute some of the charge it contains in the pres- 
ence of the field produced by the second nuclei. As long 
as distances between the bodies are large enough, we can 
disregard the strong interaction between them, and con- 
sider only the coulomb interaction. It is important to 
notice that in this context the strong force plays a role 
within each body, making the transport of charge within 
the nuclei possible. Even though the dynamics of our 
model is developed in a quantum mechanical context, it 
does not consider spin effects, nor other non classical ef- 
fects apart from those described by the scalar Schodinger 
equation. 

Our main result is that the redistribution of charge in- 
side each ball lowers the Coulomb barrier while ions ap- 
proach each other. The modification is expressed through 
the corrected differential cross sections and the Gamow 
factor, which is of special interest in astrophysical sce- 
narios (e.g., reaction networks in supernova Q)- 



II. THE MODEL 

Charge rearrangement in the presence of an external 
electric potential is common in nature. For example, 
charge move within conductors in response to external 
electric fields, as described by distance dependent image 
charges. 

In the case of nuclear reactions, there are a number 
of ways in which we can expect to obtain a charge re- 
arrangement within a nuclei, in the presence of the field 
produced by a second nuclei. For example we can imagine 
(1) a collective excitation of the nuclei, that is a fluctu- 
ation of nuclear matter around its equilibrium state; (2) 
a fluctuation involving only the surface, with an unper- 
turbed bulk, of a proton rich nuclei; (3) a local interac- 
tion affecting an unpaired proton or neutron; (4) etc. In 
all of these cases the specific description will be strongly 
dependent on the system under consideration. 

For this paper, and simplicity, lets assume that this 
charge rearrangement inside the nuclei can be modeled, 
in a first approximation, as an induced electric dipole. 
More exact calculations can also be considered, such as 
higher order moments, but for simplicity will be reported 
elsewhere. We may imagine many others ways to induce 
an electric dipole moment in a nuclei, all of which have 
their own lifetime. Of interest is the case that involves 
time-reversal violation 0, but produces an effect that 
seems to be weaker than those coming from charge re- 
arrangement 0. It is important to note, that the re- 
arrangement of charge will depend on the position and 
charge of the second nuclei. 

For the problem at hand, we can write the electric 
dipole operator as 



(1) 
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where e is the unit charge, A is the mass number, Z is 
the atomic number, and N the neutron number. The op- 
erator D has a well defined expectation value during the 
lifetime of the charge rearrangement, p = (D). Since the 
strength of any multipolar field is proportional to the size 
of the polarized object |lOl |. we will consider interactions 
of type A + p, where A is a large ion and p is a proton. 
We will neglect a possible dipole moment for the proton. 
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Suppose now that the incoming proton induce a charge 
rearrangement in the ion, observed externaUy as a dipole 
moment, which can be visuaUzed as a pair of charges 
±ge displaced a distance d away from the ion center. 
The orientation of the dipole is induced in such a way 
as to reduce the repulsion interaction, for a given sepa- 
ration r between the ion and the proton. Therefore, in 
the lhnd^o2qed = Pz, we can approximate the potential 
energy as 



U{r) 



^1^2 6^ PzZ2e 
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where pz ~ |p|, and the minus sign comes from the dis- 
cussion above. We are interested in the potential close to 
the origin, and in the limit d <C r. Since, the rearrange- 
ment of charge, as described by Pz, will depend on the 
relative position of the second nuclei, we are interested in 
situations for which pz is finite as r — > 0. This definition 
is consistent with our above assumptions. 

Due to the symmetry of the problem we choose 
parabolic coordinates: ^ — r — z, rj = r + z, 
Lp = arctan(?//x) pjj. and seek a solution of the form 



(3) 



where the function corresponds to a spherical outgoing 
wave. Operating on the scalar Schrodinger equation, we 
obtain 
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where k is the wave number, /i the reduced mass, c is the 
speed of light, v is the speed of the projectile, a is the fine 
structure constant, and is the azimuthal outgoing angle 
for the projectile (this must be small to satisfy Pz « |p|). 
This angle has meaning only for an elastic interaction. 

This equation has two linearly independent solutions 
that are regular at the origin, one for each possible sign 
of /3, so the most general solution is a linear combination 
of both 



<p{^)^C+C'^iFi{-iK + i[3,l + 2il3,ikS,) 

+ C-r'^iFi{-iK-if3,l-2ip,ik^) , 



(5) 



where iFi is the hypergeometric confluent function, and 
C± are the two constants of integration. Now we want 
to calculate the Gamow factor, which is by definition 
the quotient between the probabilities of finding close 
together particles with the same and opposite sign [T]| . 
In order to obtain a value for the constants, we normalize 
the flux asymptotically far and make C+ = C_ for sim- 
plicity. In this way we can define the following function 
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■dp — arg 
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where the function arg is the argument of the complex 
number (z = re*^, arg z = 9). From this expression we 
can write the corrected Gamow factor 
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where k± ^ and \k±\ ^ 1. Note that change from 
(3 ^ i(3 due to the change in the orientation of the dipole. 
To check this expression we take the limit (3^0, corre- 
sponding to no rearrangement of charges, and we get 



lim Gf = e"2^' 



(7) 



that is, we recover the usual expression. In this way we 
have calculated a corrected Gamow factor due to charge 
rearrangement within the target nuclei. 



III. ANALYSIS AND RESULTS 

With the aid of this simple model we can estimate in 
a first approximation how the rearrangement of charge 
in ions may modify the fusion process. As an example 
consider Fig. ^ in which we plot the Gamow window 
that is of interest in astrophysics processes ^3 • Here we 
make an estimation for a reaction IJ^^ Hi+p. In particular, 
this heavy ion has an El transition between levels with 
.r-. 7/2- and 9/2+, energy E^=214.3(keV), and half 
life 1.86(ns) jl^]. We take the ion in this excited state, 
and in order to make a simple and rough estimation, we 
consider a semi-classical calculation starting from eq.((T|l. 
which can be re-casted in the form 



NZ , 
e — — d 
A 



(8) 



where d is the distance between the center of mass for 
proton and neutron in the heavy ion. Now we take the 
energy associated to this separation d equal to E-^ as 
the work done in a Yukawa type potential This 
estimation give us Pz — 0.122(e-fm), which is not too 
far from experimental measurements in similar ions |l5j| . 
Now we need to estimate the azimuthal outgoing angle 
9, and in the same spirit, we use the Rutherford relation 
between impact parameter s and 9 
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where E is the energy of the incident particle. Taking 
E = kT, and for s the diameter of the proton in order to 
maintain 9 small, we have that for kT = 3(keV) we obtain 
9 = 0.00012(rad), while for kT = lO(keV) we obtain 9 = 
0.00039(rad). As can be seen in the figure, the impact of 
this factor is strongly dependent not only on the specific 
system considered but also on thermodynamic variables. 

Reaction cross section, which includes Gamow factor 
through the astrophysical factor, is proportional to the 
amplitude of this curve. Therefore, the rearrangement of 
charge could be important in reaction networks in stellar 
interiors, permitting the burning process of heavy ele- 
ments to start with less energy than in the case in which 
there is no charge rearrangement. 
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FIG. 1: Gamow window: The dashed line corresponds to 
the non corrected case and the soUd line represents the cor- 
rected case for = 0.122(e-fm) and 9 = 0.00012(rad). We 
use in this example the interaction ^'^^Hf+p, for kT — 3(keV). 
At the bottom the same reaction but for 9 — 0.00039(rad) and 

kT — lO(keV). (f= gf{E,/3) e~'^ , where gj is G/ as a func- 
tion of E instead of k, and have been scaled by 10^* and lO''^ 
respectively) . 



mechanism here proposed only affects the pre-compound 
nucleus stage, changing the number of such systems for 
the same incident energy but no its characteristics, leav- 
ing any subsequent nuclear evolution qualitatively unal- 
tered. 

Within the present model, the charge fraction that 
participates in the charge rearrangement, and hence the 
dipole moment Pz, can not be derived and we need to 
estimate it from other considerations (parameters Pz and 
9 given in the figures should be considered for example 
purposes only). We are developing a technique similar to 
molecular dynamic simulations (nucleon clusters) which 
may permit us not only estimate the charge fraction that 
participates in the charge rearrangement, but also it may 
provide us with the shape of the observed heavy ions as 
explained through charge rearrangement. The inclusion 
of the spin effects, which have not been considered in this 
analysis, may become relevant when taking into account 
selection rules due to angular momentum conservation. 
These considerations will be presented elsewhere. 

Now, a comment on the approximation « |p| is nec- 
essary. In eq. the position operators are subject 
to fluctuations jl7| as reflected in the vector dipole mo- 
ment, therefore, the dipole vector is not always pointing 
in the exact direction of the incoming projectile. This 
allows for the possibility of a finite A, implying an in- 
crease in the Gamow factor. On the other hand, it is 
possible that there are other contributions to the dipole 
moment, besides the induction effect. For example the 
nuclei could absorb electromagnetic radiation before the 
encounter with the proton, generating a dipole vector 
that will tend to align to the field of the incoming par- 
ticle, but not exactly because of the fluctuations. Note 
that the requirement of a small 9 ^ is compatible with 
the assumptions made in the calculation of the scattering 
cross section for the Coulomb field [Tlj . 

In the scope of nuclear astrophysics, the electric dipole 
strength in the energy region close to the neutron thresh- 
old has shown a strong influence on various processes 
of nucleosynthesis because in relevant scenarios even sub 
percent contributions to the El-energy weighted sum rule 
in this re gion may lead to completely different abundance 
patterns [I^l Of course, this has an strong connec- 
tion to p-process . 

It is not too far fetched that the proposed idea may also 
contribute in the longstanding puzzle observed in labora- 
tory measurements, in which the electron screening effect 
is surprisingly larger than theoretical prediction based on 
an atomic physics model • 



IV. CONCLUSION 



At this point, is important to notice that this charge re- 
arrangement effect is being considered in a general frame- 
work, and any detail in subsequent products will depend 
on the specific microscopic nuclear model used. For ex- 
ample, if we choose the Hauser-Feshbach formalism, the 



We have found that approximating the source of the 
coulomb field by a point may not always provide an accu- 
rate description in a fusion interaction, because it hides 
spatial degrees of freedom, and may be one source of 
disagreement between the theory and experimental mea- 
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surements. The rearrangement of charge in ions may 
tend to increase the reaction cross section, and could 
provide an extra ingredient to the understanding of some 
anomahes in reactions near the coulomb barrier without 
modification of the nuclear potential. 
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